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Abstract: Strict positive realness (SPR) is an im- 
portant concept in absolute stability theory, adaptive 
control, system identification, etc. This paper char- 
acterizes the strictly positive real (SPR) regions in 
coefficient space and presents a robust design method 
for SPR transfer functions. We first introduce the 
concepts of SPR regions and weak SPR regions and 
show that the SPR region associated with a fixed 
polynomial is unbounded, whereas the weak SPR re- 
gion is bounded. We then prove that the intersection 
of several weak SPR regions associated with different 
polynomials can not be a single point. Furthermore, 
we show how to construct a point in the SPR region 
from a point in the weak SPR region. Based on these 
theoretical development, we propose an algorithm for 
robust design of SPR transfer functions. This algo- 
rithm works well for both low order and high order 
polynomial families. Illustrative examples are pro- 
vided to show the effectiveness of this algorithm. 

Keywords: Uncertain Systems, Robustness, De- 
sign, Strict Positive Realness, Transfer Functions, 
Polynomials. 



1 Introduction 



The notion of strict positive realness (SPR) of 
transfer functions plays an important role in abso- 
lute stability theory, adaptive control and system 
identification[l-23]. Motivated by Kharitonov's sem- 
inal theorem on the robust stability for a family of 
polynomials, a number of recent papers has concen- 
trated on the strict positive realness for a family 
of transfer functions. In the spirit of Kharitonov, 
the robust SPR analysis and design problems were 
first formulated by Dasgupta and Bhagwat[8]. They 
showed that every transfer function in an interval 
transfer function family is strictly positive real if and 
only if sixteen prescibed vertex transfer functions in 
this family are strictly positive real. The sixteen crit- 
ical vertex transfer functions can be constructed ex- 



plicitly using Kharitonov's four vertex polynomials. 
This result was subsequently improved by Chapel- 
lat and Bhattacharyya, Wang and Huang, where 
only eight out of the sixteen critical vertex trans- 
fer functions need to be checked [11- 12]. For a fam- 
ily of transfer functions with affine linearly corre- 
lated perturbations, or more generally, multilinearly 
correlated perturbations, Dasgupta, Anderson et al. 
showed that it suffices to check all vertices in order 
to ensure the strict positive realness of the entire 
family [14]. By resort to the concept of positive poly- 
nomial pairs and root interlacing properties, Hollot 
and Huang solved the robust SPR design problem for 
low order and structured families [9- 10]. Anderson et 
al. considered the general robust SPR design prob- 
lem, and by using the Hilbert transform, provided 
a constructive method[14]. Betser and Zeheb made 
some further improvements [15]. 

This paper characterizes SPR regions in coefficient 
space and presents a robust design method for SPR 
transfer functions. We first introduce the concepts of 
SPR regions and weak SPR regions and give a com- 
plete characterization of them. We show that the 
SPR region associated with a fixed polynomial is un- 
bounded, whereas the weak SPR region is bounded. 
We then prove that the intersection of several weak 
SPR regions associated with different polynomials 
can not be a single point. Furthermore, we show 
how to construct a point in the SPR region from a 
point in the weak SPR region. Based on these the- 
oretical development, we propose an algorithm for 
robust design of SPR transfer functions. This algo- 
rithm works well for both low order and high order 
polynomial families. Illustrative examples are pro- 
vided to show the effectiveness of this algorithm. 



2 Preliminaries 

Denote P n as the ra-th order real polynomial family, 
R n as the n dimensional real field, and H n C P n as 



the set of all n-th order Hurwitz stable polynomials. 

In the following definitions, b(-) G P m , a(-) G P n , and 
p(s) = b(s)/a(s) is a rational function. 

Definition 1 p(s) is said to be strictly positive 
real(SPR),denoted as p(s) GSPR,if b(s) G P n ,a(s) G 
H n , and Re[p(jw)] > 0, Vw e R. 

Definition 2 p(s) is said to be weak SPR 
(WSPR), denoted as p(s) eWSPR,if b(s) G 
P"-\a(s) e H n , and Re[p(jw)] > 0, Vw e R. 

Definition 3 Given a(s) G H n , the set of the co- 
efficients (in R n+1 ) of all the 6(s)'s in P n such that 



ffl ( s ) = s " + fllS "- 1 + -.. + fflll eff' 1 (1) 

Denote as Sli a the set of the coefficients of all 
the b{s) = s n + x^- 1 + ■■■ + x n G P n , i.e., 



(xi, X2, ■ ■ ■ , x n ) in R n , such that p(s) 



6(f) 
a(s) 



GSPR; 



and denote as Sl]^ the set of the coefficients of all 

the b(s) = S™- 1 + x lS n - 2 + ■■■ + x n - X G P n_1 , 

i.e., (xi, ^2, • • • , £n-i) in , such that p(s) = 

b(s) 



a(s) 



GWSPR. Obviously, we have 



p(s) 



a(s) 



GSPR is said to be the SPR region {1} x O lo = ( (1, — , — , ■ 

I ' fro ' fro ' 



associated with a(s), denoted as tt a . 

Definition 4 Given a(s) G i? n , the set of the co- 
efficients (in R n ) of all the b(s)'s in P™ _1 such that 



bn 
fro 



p(s) := 



Kf) 
a(«) 



GWSPR is said to be the WSPR region 



associated with a(s), denoted as Sl^ 7 . 

For notational convenience, Q, a (Q^) sometimes also 
stands for the set of all the polynomials b(s) in 

pnrpn-l^ guch ^ . = Kf) G SPR(WSPR). 

o(s) 

From the definitions above, it is easy to get the fol- 
lowing properties: 

Property If p ( a ) eSPR(WSPR), then 

7T 

|cw"5(fr(iw))-arc/(a(iw))| < - , G i?, where arg(-) 
stands for the argument of the complex number, and 
the difference of two arguments can differ by an in- 
teger number of 2tt. 

Property 2^ 9 ^ Given a(s) G H n , ft a is a non- 
empty, open, convex cone in R n+1 . 

Property 3I 10 ' 11 ] Given a(s) G H n , we have fl a C 

H n ,nf c h ,i -\ 

The problem we are interested in is: Given a fam- 
ily of Hurwitz stable polynomials, how can find a 
fixed polynomial such that their ratios will be SPR- 
invariant? In what follows, we will first give some 
characterization of WSPR regions, and then propose 
an efficient design procedure for this problem. 



V(& ,fri,fr 2 ,--- A) efi a (2) 



fri b 2 6 n _i, 



fro ' fro ' ' fro 

V(&o,&l,&2,"-,&n-l)eft^ ( 3 ) 



For notational convenience, Qi a (Sl]^) sometimes also 
stands for the corresponding polynomial set. 

As we know! 9 ' 10 !, Sl a is a non-empty, open, convex 
cone in R n+1 . Thus, £l a is an unbounded set in 
R n+1 . In what follows, we will show that Sli a is 
also an unbounded set in R n . 

Theorem 1 Given a(s) G H n , ili a is a non-empty, 
open, unbounded convex set in R n . 

Proof Obviously, we have a(s) G Sl a . If the leading 
coefficient of a(s) is do, then 



ai a 2 a„ 



Thus, Sli a is not empty. 

Moreover, (l,x 1 ,x 2 , ■ ■ ■ ,x n ) G Sl a , V(xi, x 2 , ■••,£«) G 
Sli a . By Property 2, £! a is open. Thus, 
these exists 5 > 0, such that, when 
- yo) 2 + (xi - yi) 2 + ■■■ + (x n - y n ) 2 < S, 
we have (y , y 1 , y 2 , ■ ■ • , y n ) G ^a- For this <5, if 
(zi,Z2, ••• , 2r») £ P™ satisfies 



3 Geometric Characterization of 
SPR Regions 

By definition, an SPR (WSPR) transfer function 
times a positive integer is still SPR (WSPR). Thus, 
without loss of generality, let 



sj (Xl - Zl) 2 H h (x n - Z n ) 2 < S. 

Then, we have 

y/(l - l) 2 + (an - z x ) 2 + ■ ■ ■ + (x n - z n ) 2 < 5. 

Hence (1, z\, z 2 , ■ • • , z n ) G Sl a . Thus, we have 
(zi,Z2, • • • , z n ) G Sli a . Namely, Q,\ a is an open set. 



By definition, la is convex. 

In what follows, we will prove that £l la is unbounded. 
For this purpose, we first introduce some notations, 
which are needed in other proofs as well. Let 

a(s) = s n + a 1 s n - 1 +••• + <*„ G H n (4) 



b(s) = x s n + x^ 1 - 1 + --- + x n eP n U P n - X (5) 
Then Vw £ R, we have 



Re[ 



^4] = T7 ^MHMa(-M 

a(ju>) Hjuj)\ 2 



i«-)P ti 



where q := ^ a-kX2i-k{— where ao = 1, and 
let di = 0, .Ti = 0, when i < or i > n, I = 0, 1, • • • , n. 



Introducing the matrices 





ai 


1 








• 


• " 




03 


a 2 


ai 


1 


• 


• 




a 5 




«3 


a 2 


ai • 


• 






1 CL2n- 


2 «2n- 


3 &2n- 


4 a 2 n-5 ' 


• a n _ 



En '■- 



A :-- 



-1 



1 

-a 2 

0,4 E n H a E n 



X 
X\ 

Xn—1 



b :-- 



c := 



co 

Cl 



where ai 
that 



when i > n. Then, it is easy to verify 



c = Ab 



(6) 



Since a(s) € H n , we know that E n H a E n is in- 
vertible. For any d — [di, d,2 1 • • ■ , d n ] T G R n 
such that all elements of d are positive. De- 
note a = [— a2, &4, — a6, as, • • • , (— l)"a2n] T , where 
ai = when i > n. Let b = (E n H a E n )^ 1 (d - 
a) :— [bi, b 2l ■ ■ ■ , b n ] T , Then obviously, we have 
[l,d 1 ,d 2 , ■ ■ ■ ,d n ] T = A[l,b 1 ,b 2 ,---,b n ] T , and b n = 

— . By c = Ab, we know that (b\,b 2 , • • • , b n ) e Qi a . 

On the other hand, due to the arbitrariness of d, 
d n can be taken arbitrarily large. Therefore, b n can 
also be arbitrarily large. Namely, Qi a is unbounded. 
This completes the proof. 

Remark 1 Given two stable polynomials ai(s) and 
02(5) , existence of a polynomial b(s) such that 
-^A and are both SPR is tantamount to non- 

ai(s) a 2 (s) 

emptiness of the intersection of the two SPR regions 
associated with &i(s) and 02(5) . Since f2 a and f2 la 
are both unbounded sets. When dealing with robust 
SPR design problem, we must find the intersection 
of several unbounded sets (i.e., SPR regions), which 
is intractable. This is the reason that we introduce 
the concept of WSPR regions, which are bounded as 
shown below. 

Theorem 2 Given a(s) £ H n , 0]^ is a non-empty, 
bounded convex set in R n . 

Proof 

H a , E n and A were defined in the proof of Theorem 
1. 

Denote B as the (n — 1) x (n — 1) matrix formed by 
the first n — 1 row and last n — 1 column of the matrix 
E n H a E n . Obviously, £> is also invertible. 

Denote a := [a\, —03, 05, —07, 
( ai = when i > n). 

Let 6 := — i? _1 a = [61, 62, • • • , 6„_i] T . Since a(s) € 
it is easy to verify that & n _i > 0. Denote b = 
[0, 1, bi, 6 2 , • • • , &n-i] T - Let c = ci = • • • = c n _i = 0, 
c n = a n b n -i, in c := [c , Ci, • • • , c n ]. Then, it is 
easy to verify that c = Ab is true. Thus, we have 
(61, 6 2 , • • • , &„_i) G Namely, is not empty. 

By definition, is convex. 

We now prove that is bounded. 



(-l)"- 1 a 2(n _ 1)+1 ] r 



For any x 2 , • • • , £„-i) E 0]^, we have 



la' 

s"- 1 +xi5"- 2 + --- + a: K _ 1 
s" + ais"- 1 + --- + a„ 

Moreover, Vw € -ff, we have 



e WSPi?, 



By Property 3, s"" 1 + xis"" 2 + • • • + a; n _i £ H 



n-l 



Re( 



Thus 



Xn—1 



n-l 



a,. 



=*-) > 0. (7) 



5" + ais"" 1 

^ c n— 1 _i_ . 7 i H c n— 2 



— 1 



8 =j W ) > 0, e R. 

(8) 



Obviously 



■ ais"^ 1 + ■ • 



= s + 



s n — 1 + xis n 2 + ■ • • + 

(ax - ii)s" _1 + (a 2 — x 2 )s" -2 + • ■ • + (o„_i - x„_i)s + a„ 



s" 1 + lis" 2 + ■ ■ ■ + cc n _i 



Therefore 



Re{ 



(oi — xi)s n 1 + (a 2 — x 2 )s" 2 + ■ • ■ + (o n _i — x n _i)s + a„ 
s 71-1 + XiS" -2 + • ■ ■ + X n -1 



_ , »" + »is" H h On i > „ , . n 

s n + ais" -1 + ■ • • + a n 

:Re (^TT ^rj- t U=j u )>0, Vc^efl. 

s n 1 + ^is n + - • • + x n — i 



C {(a;i,a;2,---)^n-i)|0 < x 1 < a 1: < x 2 <a 2 ,---, 

< x n -! < a n _i}. 

Thus, is bounded. This completes the proof. 

Remark 2 It should be pointed out that is not 
an open set in i? n_1 . In fact, from the proof of The- 
orem 2, we know that is tangent to the hyper- 
plane xi = a\ in R 71-1 . And there exist some points 
of VlY a in this hyperplane. Thus, Qj^ can not be an 
open set. Obviouslyf^ is a non-empty, convex cone 
in J? n_1 . Thus, is also unbounded. One may be 
tempted to believe that is not an open set either. 

Though VtY a is not an open set. The following theo- 
rem guarantees such a fact: when the intersection of 
two or more WSPR regions is not empty, then the in- 
tersection must be a region, not a single point. This 
means that Ackermann's counterexample (that the 
unstable region is an isolated point so that gridding 
the parameter space can lead to erroneous conclu- 
sions no matter how dense the gridding is[19]) does 
not happen in this case. 

t) Theorem 3 Given a(s) E H n , if 
(xi, X2, ■ ■ ■ , x n -\) E O]^, then for sufficiently small 
e > 0, we have (xi — e,x 2 — e, - ■ ■ , x n -i — e) E f2^. 

Proof V(xi, x 2 , ••• , x n -\) E SlYai an( i ^ w e we 
have 



It is easy to see that 



(ai — xi)s n 1 + (a 2 - X2)s" 2 + • ■ • + (o n -l — x„-i)s + a„ 
a"- 1 + us"" 2 H h z„-i _ ~ 



e {SPR} U {WSPR}. 



Again, by Property 3, we have 
(01 - xxjs 11 - 1 + (a 2 - x 2 )s n - 2 + 



+ (a n -i-x n -i)s + a n eH n - 1 UH n - 2 (9) 



Hence 



< x\ < a\, < x 2 < a 2 , ■ ■ ■ , < x n -\ < a„_i 

(10) 

Namely 



QZ C {{x llX2l - ■ ■ ,x n -i)\a{s) := 

n 

^(a 4 - a*)*" - ' E H n ~ l U iT l ~ 2 , where jc n = 0} 
»=i 



Re( - ; - — — j — |»=ja,) > 0. (11) 



s™ + ais™- 1 + • • • + a n 
Ve > 0, since 



Re{- 



n-l 



+ (xi - e)s"- 2 + • • • + (s„_i - e) 



s n + ais 



n-l 



H ha„ 



s™ + 01S™- 1 H ha r , 



+i?e( 



(-£)(* 



n-2 



c.n-3 



+ ••• + 1) 



= Re( 



s n 1 + a;is 



s™ + ais"" 1 

n-2 i 



H hfl„ 

^n-1 , x 



\s=ju>) 



Or, 



+ ■ 



a{juj)\'< 



(V^ + ^H), 



where c(u>) is a real polynomial of order less or 
equal to 2 (n — 2). Thus, when | uj \ is sufficiently 
large, the sign of (— e)(— cj 2 ^ n ^ 1 '> + c(w)) will be pos- 
itive. Namely, there exists u)\ > such that, for all 

uj |> 0J 1 , 



Obviously, b(s) and a(s) have same order n. Thus, 
s n-i _|_ r Xl _ e \ s n-2 _|_ |_ r x _ i _ £ \ there exists ui 2 > such that, for all | w |> W2, we 

^ a n +ai „n-l " !-=*-)> 0. 

s +ais H ha„ have Re( 7 , ) > 0. 

(12) a(jw) 

Denote Denote 



Mi= inf Re( — + ^ . + ' - + | s= ,^), 

(13) 



Nx = sup | Re( 



|u,|<u,! \a(ju 



{uj 2{n -^ - c(w)) | (14) 



Mi 

Then Mi > and iVi > 0. Choosing < e < — , 

iVi 

by simple computation, we have 

Re( s "' 1 + (3 ; 1 7 £)s I'; + --- + (3: "- 1 - £) U^) > o.Vu, e R 

s n + ais n 1 + ■ ■ ■ + a n 

(15) 

Therefore, 

^- i +(^-^-;+-+(^-i- e ) e ws 

s" + ais™- 1 H h a„ 

namely (x\ — e,X2 — £,■•■, Xn-i — e) € f2^. This 
completes the proof. 

The following theorem shows the relationship be- 
tween SlYa an d ^a, and plays an important role in 
robust SPR design. 

Theorem 4 Given a(s) G H n , if 
(zi,a;2, • • - ,a; n _i) £ Q}£, then V(l,ai,a 2 , ■ ■ ■ ,a n ) G 
i? n+1 , we can take sufficiently small e > such that 
(0, l,x 1 ,x 2 ,-- -,x n -i) +e(l,ai,«2, ■■■ ,a n ) E o . 
Proof Denote 



M 2 = inf Rc(^l), N 2 = sup | Re(^M) | 

(21) 
M 2 

Then M 2 > and N 2 > 0. Choosing < e < — -, 

jv 2 

by simple computation, we have 



Re(^4) > 0,Vco G R 



(22) 



This completes the proof. 



4 Applications in Robust Design of 
SPR Transfer Functions 

Generally speaking, the design problem is more dif- 
ficult than analysis problem, since it is usually con- 
structive, i.e., it not only shows the existence of the 
solution, but also provides a constructive procedure 
to find it. In this section, we will propose an algo- 
rithm for robust design of SPR transfer functions. 
This algorithm works well for both low order and 
high order polynomial families. Illustrative exam- 
ples are provided to show the effectiveness of this 
algorithm. 

Suppose 



b( S ) = S™- 1 + X lS n - 2 + ■ ■ ■ + (16) 

a(s) = s n + a 1 s n - 1 + --- + a n , (17) 

b{s) = b(s)+ea(s). (18) 
Since (x\, x 2 , • • • , £ &iai we have 



Re (MM) >0 ,Vu,Gi*. 



(19) 



We only need to show that, for sufficiently small 
e > 0, 



Re(^4) > 0,Vco G R 



(20) 



{a,i(s) 



{i) s n -\i= l,2,---,m.} 



(23) 

How can we find a polynomial b(s), such that Pi(s) := 
eSPR,i = l,2,---,m? 

ai(s) 

As observed earlier, existence of such a polynomial 
b(s) boils down to the condition that the intersection 
of the SPR regions associated with <Zj (s) is not empty. 
From the results in the previous section, we know 
that SPR regions are unbounded, whereas WSPR re- 
gions are bounded. Thus, by a computational consid- 
eration, we first consider the intersection of WSPR 
regions, and then construct a polynomial b(s) by us- 
ing the technique presented in the previous section. 



Since SPR(WSPR) transfer functions with fixed nu- 
merator (or denominator) enjoy convexity property, 
namely, if there exists a polynomial c(s), such that 

and ^4 are b °th SPR(WSPR),thcn, it is easy 
a(s) b(sj 

to verify that, for any a > 0, (3 > and (a, f3) ^ 
c(s) 

(0, 0), we have — . . mi \ eSPR(WSPR). There- 
y ' aa(s) + f3b(s) y ' 

fore, the assumptions made on F do not lose any gen- 
erality. Actually, the method proposed in our paper 
also applies to convex combination of polynomials, 
interval polynomials, and more generally, polytopic 
polynomials and multilinearly perturbed polynomi- 
als[l,4,14]. 

By the results presented in the previous section, we 
propose the following design procedure: 

Step 1. Test the robust stability of the convex hull 
of F, i.e., F. If F is robustly stable, then go to Step 
2; otherwise, print "there does not exist such a b(s) 
" ; (by Definitions 1 and 2) 

Step 2. Let a; = min{a[ l \i = 1, 2, ■ • • , m}, I = 
1, 2, • • ■ , n — 1. Gridding the hyperrectangle 

D:= 

{{x 1 ,x 2 ,- ■ ■ ,x„_i) | < xi < a u l = 1,2,-- -,n- 1} 

(24) 

according to the precision required; (by Theorem 2 
and its proof) 

Step 3. Take b := (&i, 6 2 , • • • , & n -i) at eacn gridding 
point. Test if b belongs to n^O^. by the following 
steps: 

1) Test if the (n — l)-th order polynomial with coeffi- 
cients (1, &i, &2) ■ " ' > b n -i) belongs to H n ^ 1 (by Prop- 
erty 3), 

2) For i — 1, 2, • • • , m, test if the polynomial with 
coefficients (a^ — bi,a 2 ^ — b 2 , • ■ ■ , — &n-i, On^) 
belongs to UH n ~ 2 , respectively (by Theorem 
2 and its proof), 

3) Test if b belongs to 0]^. , i = 1, 2, • • • , m . 

If all 1), 2), 3) above are satisfied, go to Step 4; oth- 
erwise, move to the next gridding point and test 1), 
2), 3) again; (If all gridding points have been tested, 
then print "there does not exist such a & in n™ 
with the given precision" ) . 

Step 4. Take a sufficiently small e > 
such that (£r,l,6i,&2,-",&n-i) G ^T=i n a % - 
Hence, the n-th order polynomial with coefficients 
(e, 1, 61, 6 2 , • • • , b n -i) satisfies the design require- 



ment, (by Theorem 4). 

For the low order stable interval polynomial family 
or low order stable convex combination, existence of 
the solution to the design problem is always guaran- 
teed[15,16]. Given adequate precision, our method 
will surely find a polynomial that satisfies the de- 
sign requirement. As shown by numerous examples 
below, our method is also effective for higher order 
polynomial families. 

Example 1 

Let 

F = { ai(s) = s 4 + lls 3 + 56s 2 + 88s+l, 
a 2 {s) = s 4 + lis 3 + 56s 2 + 88s + 50, 
03(5) = s 4 + 89s 3 + 56s 2 + 88s + 1, 
a 4 (s) = s 4 + 89s 3 + 56s 2 + 88s + 50} 

(25) 

the methods proposed in [8-10,14,15] do not work 
here. Using our method, it is easy to get b(s) = 
s 3 + 3.3s 2 + 2.24s + 1.76 £ n 4 =1 0^. Then let 
c(s) := es 4 + b(s), where e > is sufficiently small, 
e.g., let e < 0.3 (which is determined by Theorem 
4), it is easy to check that the design requirement has 
been met. 

Note that the example above is constructed by over- 
bounding the line segment in [13] by an interval poly- 
nomial family. Thus, instead of dealing with two 
vertex polynomials as in [13], we must now deal with 
four Kharitonov's vertex polynomials. 

In what follows, we will give two more examples of 
higher order polynomial families. 

Example 2 

Let 

F = { 

01 (s) = s 7 + 9s 6 + 31s 5 + 71.5s 4 + 111.5s 3 + 109s 2 + 76s + 12.5, 
a 2 (s) = s 7 + 9.5s 6 + 31s 5 + 71s 4 + 111.5s 3 + 109.5s 2 + 76s + 12, 
a 3 (s) = s 7 + 9s 6 + 31.5s 5 + 71.5s 4 + Ills 3 + 109s 2 + 76.5s + 12.5, 
a 4 (s) = s 7 + 9.5s 6 + 31.5s 5 + 71s 4 + Ills 3 + 109.5s 2 + 76.5s + 12} (26) 

It is easy to see that the convex hull F of F is 
robust stable. Using our method, it is easy to get 
b(s) = s 6 + 7.2s 5 + 18.6s 4 + 42.6s 3 + 44.4s 2 + 43.6s + 
15.2£n 4 =1 ^. 

Then let c(s) := es 7 + b(s), where e > is sufficiently 
small, e.g., let e < 0.1, it is easy to check that the 
design requirement has been met. 

Example 3 



Let 

F = { oi(s) = s 9 + lis 8 + 52s 7 + 145s 6 + 266s 5 + 331s 4 

+280s 3 + 155s 2 + 49s + 6, 
a 2 (s) = s 9 + lis 8 + 52s 7 + 146s 6 + 265.5s 5 + 332s 4 

+278. 5s 3 + 151s 2 + 48s + 2} 
(27) 

It can be verified that a 2 (s) — a\(s) = s 6 — 0.5s 5 + s 4 — 
1.5s 3 — 4s 2 — s — 4 satisfies the extended Alternating 
Hurwitz Minor Condition[ 7,21,22 ]namely, it is a con- 
vex direction for Hurwitz stability! 7,21,22 ]. Moreover, 
it is easy to see that a\ (s) and a 2 {s) are both Hurwitz 
stable polynomials. Thus, the convex hull F of F is 
robust stable' 7,21,22 '. Using our method, it is easy 
to get b(s) = s 8 + 8.8s 7 + 41.6s 6 + 87s 5 + 159.3s 4 + 
132.4s 3 + 111.4s 2 + 30.2s + 9.6 G Qf ai n ^Z 2 ■ Thus, 
let c(s) := es 9 + b(s), e > 0, e sufficiently small, e.g., 
take e < 0.07, then the design requirement has been 
met. 

Note that our design method is also effective when 
dealing with discrete time systems. Note also that, 
in the Examples 1-3, b(s) is not unique. Using our 
method, we can get all such 6(s)'s with given preci- 
sion. 

It should be pointed out that there is hardly any ex- 
ample with order higher than 4 in the literature. Re- 
cently, a sixth-order example of interval family was 
given in [23] as follows. Unfortunately, this example 
is incorrect. 

Example 4 

Suppose 

f = { 

oi(s) = s 6 + 0.8s' 5 + 58.06s 4 + 50.9s 3 + 1028.5s 2 + 163.82s + 1042.5, 
a 2 (s) = s 6 + 1.5s 5 + 58.06s 4 + 28.3s 3 + 1028.5s 2 + 376.36s + 1042.5, 
a 3 (s) = s 6 + 0.8s 5 + 68.62s 4 + 50.9s 3 + 755.47s 2 + 163.82s + 3286.7, 
a 4 (s) = s 6 + 1.5s 5 + 68.62s 4 + 28.3s 3 + 755.47s 2 + 376.36s + 3286.' 

(28) 

find a polynomial b(s), such that 

Pi (s) :=-^Sr £SPR,i = 1,2,3,4. 

By Definition 1 , Definition 2 and Property 3, a pre- 
requisite of the robust SPR design problem is that 
the convex hull F of F is robustly stable. But 
it is easy to check that F is not robustly stable. 
In fact, (1.0446 ± 5.8969i) are roots of oi(s) with 
positive real part; (1.037 ± 4.9835i) are roots of 
a 2 (s) with positive real part; (0.03291 ± 7. 5026i) and 
(0.68089 ± 2.4933i) are roots of a 3 (s) with positive 
real part; (0.87123 ± 2.867i) are roots of a4(s) with 



positive real part. Thus, it does not make sense to 
consider the robust SPR design in this case. 

It should also be pointed out that, for the ver- 
tex set of a general polytopic polynomial family 
F = Ms) = s" + Er=i°l (<) = l,2,---,m.}, 

even if F is robustly stablest is still possible that 
there does not exist a polynomial c(s), such that, 
c(s)//(s) eWSPR, for all f(s) e F. Namely 

n™ ^ = 0. 

To see this, let us look at an example of a third order 
triangle polynomial family. 

Proposition ll 15,16,18 ! 

Let a(s) — s 3 + ais 2 + a 2 s + a 3 6 H 3 , then 

= {(x, y)\ajx 2 + 2(2a 3 - ai a 2 )xy + a\y 2 
—2a 2 a 3 x — 2aia 3 y + a 2 < 0} 
U{(x,y)\x < a-L,a 2 x - a x y - a 3 > 0,y > 0}. 

Example 5 

Let 



F = { ai(s) = s 3 + 2.6s 2 + 37s + 64, 

a 2 (s) = s 3 + 17s 2 + 83s + 978, (29) 
a 3 (s) = s 3 + 15s 2 + 28s + 415} 

it is easy to verify that ai(s),i — 1,2,3, are Hur- 
witz stable. Moreover, all edges of F, i.e., Aa^(s) + 
(1 — A)aj(s),A G [0,1], i,j = 1,2,3, are also Hur- 
witz stable. Therefore, by Edge Theorem! 1,6,7,20 ], 
F is robustly stable. On the other hand, by a 
direct computation using Proposition 1, we have 
nflZ nnZ = 6. Henceforth, there does not ex- 

Ld\ 1(22 LG.3 ' ' 

ist a polynomial c(s) such that c{s)/ai{s) 6 WSPR, 

t = i, 2, 3 (although nf ai nnZ 2 + K ai + 4>, 
and n% 2 nnZ 3 ^<f>)- 

Note that, in this example, though we have (~l 
Q]^ 2 (~1 ^Ya 3 — 4>i we do not know whether fi 0l (~1 
O a2 n O a3 is an empty set or not. This is a problem 
deserving further study. 

For a fourth ( or higher ) order stable interval poly- 
nomial family ( or stable convex combination of two 
polynomials), does there exist a polynomial such that 
their ratios are SPR-invariant? This is still an open 
problem! 1,8,9,13,14 ]. From our numerous examples, it 
seems that such a polynomial can always be found. 



Thus, we conjecture that this problem has a positive 
answer. 
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